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We discuss a hydrodynamical description of the eigenvalues of the Polyakov line at large but finite 
Nc for Yang-Mills theory in even and odd space-time dimensions. The hydro-static solutions for 
the eigenvalue densities are shown to interpolate between a uniform distribution in the confined 
phase and a localized distribution in the de-confined phase. The resulting critical temperatures 
are in overall agreement with those measured on the lattice over a broad range of Nc, and are 
consistent with the string model results at Nc = oo. The stochastic relaxation of the eigenvalues 
of the Polyakov line out of equilibrium is captured by a hydrodynamical instanton. An estimate of 
the probability of formation of a Z(Nc) bubble using a piece-wise sound wave is suggested. 

PACS numbers: 12.38Aw, 12.38Mh, 71.10Pm 


1. Introduction. Lattice simulations of Yang-Mills 
theory in even and odd dimensions show that the con¬ 
fined phase is center symmetric mm- At high temper¬ 
ature Yang-Mills theory is in a deconfined phase with 
broken center symmetry. The transition from a center 
symmetric to a center broken phase is non-perturbative 
and is the topic of intense numerical and effective model 
calculations [3] (and references therein). Of particular in¬ 
terest are the semi-classical descriptions and matrix mod¬ 
els. 

In the semi-classical approximations, the confinement- 
deconfinement transition is understood as the breaking of 
instantons into a dense plasma of dyons in the confined 
phase and their re-assembly into instanton molecules in 
the deconfined phase n m- A mechanism similar to 
the Berezinsky-Kosterlitz-Thouless transition in lower di¬ 
mensions [6], and to the transition from insulators to 
superconductors in topological materials [7]. In matrix 
models, the Yang-Mills theory is simplified to the eigen¬ 
values of the Polyakov line and an effective potential is 
used with parameters fitted to the bulk pressure to study 
such a transition Eiin], in the spirit of the strong coupling 
transition in the Gross-Witten model m- 

Matrix models for the Polyakov line share much in 
common with unitary matrix models in the general con¬ 
text of random matrix theory HU. The canonical exam¬ 
ple is Dyson circular unitary ensemble and its analysis 
in terms of orthogonal polynomials or a one-component 
Coulomb plasma. The Dyson circular unitary ensemble 
relates to the Calogero-Sutherland model which is an ef¬ 
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fective model for quantum Luttinger liquids. A useful 
analysis of this model uses the collective quantization 
method developed in m with its hydrodynamical inter¬ 
pretation [Hill]. 

In this letter we develop a hydrodynamical description 
of the gauge invariant eigenvalues of the Polyakov line 
for an SU(Nc) Yang-Mills theory at large but finite Nc- 
We will use it to derive the following new results: 1/ a 
hydrostatic solution for the eigenvalue density that inter¬ 
polates between a confining (uniform) and de-confining 
(localized) phase; 2/ explicit critical temperatures for the 
Yang-Mills transitions in 1 -|- 2 and 1-1-3 dimensions; 3/ 
a hydrodynamical instanton for the density distribution 
that captures the stochastic relaxation of the eigenval¬ 
ues of the Polyakov line; 4/ an estimate of the fugacity 
or probability to form a Z(Nc) bubble using a piece-wise 
sound-wave. 

2. Polyakov line in 1 -|- 2 dimensions. The ma¬ 
trix model partition function for the eigenvalues of the 
Polyakov line for SU(A),) in 1 -|- 2 dimensions was dis¬ 
cussed in |5]. If we denote by diag(e*®D •■•i ) with 

y). dj = 0 the gauge invariant eigenvalues of the Polyakov 
line, then |S] 

Nc Nc 

Z[a,P] = ( 1 ) 

i=i i<j 

with Zij = Zi — Zj and Zi = e*®'. The perturbative po¬ 
tential V{zij) is center symmetric and quadratic in lead¬ 
ing order or V{\zij\) « with a{T) = T‘^V2/2'k 

and V 2 the spatial 2-volume |8]- The mass expansion 
of the one-loop determinant gives /3(r) = V 2 /tt |8]- 
The Debye mass is self-consistently defined as = 




2 


Ncg'^T(hi{T/mD) + C')/27r [T^ to tame all infra-red di¬ 
vergences, with C « 1.3 from lattice simulations [niiiT]. 

Q can be regarded as the normalization of the squared 
and real many-body wave-function which is the 

zero-mode solution to the Shrodinger equation i/o'I'o = 0 
with the self-adjoint squared Hamiltonian 


dtV = i[H,v\= ( 8 ) 

-de -I- - deA - Aclglnp -I- T:j5[Ap\H - 2a[Ap]s) 

with the sine-transform [Ap]s = / sin(0 — 0')A{9')p{9'). 
Note that all the relations hold for large but finite 


ATc 

Hq = 'y ^ {—di -\- a^) [di + B 4 ) (2) 

i=l 

with di = d/89i and the pure gauge potential = diS. 
Here = — lnd>o[ 0 ] is half the energy in the defining 
partition function in 0. In 0 the mass parameter is 
1 / 2 . 

3. Hydrodynamics. We can use the collective co¬ 
ordinate method in [T^] to re-write 0 in terms of the 
density of eigenvalues as a collective variable p[9) = 
~ ^i)- I^or that, we re-define Hq ^ H through 
a similarity transformation to re-absorb the diverging 2 - 
body part induced by the Vandermond contribution A = 
nz<j i-e- ^ = '^o/'/A and = Hq^AA. 

Now H is of the general form discussed in m and is 
amenable after some algebra to 


H = J d9 {deirpdeTT + pu[p]) (3) 

with the potential-like contribution 

u[p] = (^A(0) - ^ (4) 

Here 

Ai9) = ia(r) J d9'p{9') deV (^2 sin (^) ) (5) 

and ph is the periodic Hilbert transform of p 

[p\h = Ph[0) = j p( 6'') cotan (6) 

As conjugate pairs, 7 r( 0 ) and p{9) satisfy the equal-time 
commutation rule ['k{ 9), p{9')\ = —i{5{9 — 9') — l/2-K). 
We identify the collective fluid velocity with v = dgn and 
re-write 0 in the more familiar hydrodynamical form 


4. Hydro-static solution. The static hydrody¬ 
namical density follows from the minimum of 0 with 
v{9) = 0, 


I3{T)t,ph (9) - dg\np{9) = 2A{9) (9) 

To solve 0, we insert the leading quadratic contribution 
A{9) « 2d(r)sin2(6»/2) in 0 


PPh — adgp = bcipsui{9) (10) 

with a = 1 /ttP{T), b = 2a{T)/fj{T) and ci the first 
moment of the density or ttci = p{9)cos9d9. Let 

Po = A’c/27r be the uniform eigenvalue density and pi = 
p — Po its deviation. Consider the Cauchy transform 


G{z) 



Piiv) 

T] — Z 


dr] 


(11) 


with 77 = e*®. The contour C is counter-clockwise along 
the unit circle. G{z) is a holomorphic function in the 
complex z-plane. Let G~^ and G~ be its realization inside 
and outside C respectively, so that 


G^{z^e^<^) = ±Pii9) + ipH{9) 


( 12 ) 


We now carry the Hilbert transform on both sides of ( |10[ ) . 
Setting G{z) = G+{z) and using 2[piPh]h = Ph 
we have for ( 10 ) 


-b (po - ^bci{z - z ^))G + azd^G = bcipoz + ^bcf 


(13) 


on the boundary C, thus within the circle. Here, we 
should require G(z = 0) = 0 to ensure that pi integrates 
to zero. 


a « I/V 2 is subleading and will be dropped. Thus (13) 
is algebraic in G{z). Since p{9) = po + ReG+( 2 : = e^, 
careful considerations of the singularity structures of the 
quadratic solutions to (131 yield (0 is a step function) 


iL « y d9 p{9) (v'^ + u[p]) J d9p{9)\v + iAf (7) 

modulo ultra-local terms. The Heisenberg equation for p 
yields the current conservation law dfp = —2dg (pv), and 
the Heisenberg equation for v gives the Euler equation 


p{9) = \/bci{cos9 + 1) 2 {cos9 — cos9o) 2 0(|0ol ~ 1^1) 

(14) 

The analytic properties of G{z) fix ci/po = 1-|-(1 —l/^)® 
and 9q at cos^o = l~2po/5ci. For 6 < 1 the non-uniform 
solution with pi 7 b 0 is absent. For 6 ^ 1, ci —>■ 2po and 
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p{e) ^\/ 86-46202 ( 15 ) 

ZTT 

Therefore ( [T4| interpolates between a uniform density 
distribution po (confined phase) and a Wigner semi¬ 
circle (deconfined phase) with a transition at 6 = 1 or 
Tc = niD- In 1 -|- 2 dimensions the fundamental string 
tension is given to a good accuracy by = 

((1 —l/iV^)/87r)2 [5D]. Thus the ratio in 1-1-2 dimensions 


Tc C / Stt \ 



(16) 


with C « 1.3 linillT]. In Fig. [^we show the behavior 
of (16) (upper curve) versus Nc, in comparison to the 
numerical fit Td = 0.9026 -I- 0.880/fVf to the lattice 
results (lower curve) in [21]. Amusingly, ( |I^ at large Nc 
is consistent with in the string model [H]. 





The /? contribution is the self Coulomb subtraction and is 
consistent with the subtraction in the Hilbert transform. 
The saddle point equation SV j^p — ^ following from 


191 is in agreement with the hydro-static equation 


d (5r[Q!, /3; p\ 
M Sp{0) 


= 2A = 0 


( 20 ) 


6. Hydrodynamical instanton. The fixed time 
zero energy solution to 0 is an instanton with imaginary 
velocity v = —iA. We have checked that this is a solution 
to ^ for all times. The current j = pv = —ipA is 
conserved.Thus drP — 2dg{pA) = 0 or 


drP + ^{T)de{TrppH) = dgp -f 2de{pA{9)) (21) 


for Euclidean times r = it. For A = 0 and /3(T) = 2, (21) 


agrees with the viscid Burger's equation describing large 
Wilson loops in 1 -I- 1 dimensions [19] . Following M we 
identify r with the stochastic (Langevin) time. (21) de¬ 
scribes the stochastic relaxation of the eigenvalue density 
of the Polyakov line (out of equilibrium) to its asymptotic 
(in equilibrium) hydro-static solution. 


7. Sound waves. The hydrodynamical action fol¬ 
lows from standard procedure.The momentum n{0) = 
{l/de)v is canonically conjugate to the density p, and 
the Lagrange density is L = irdtp — H. Thus the action 
S = J dtdO p{9) {d — u[p]), which is linearized by 


FIG. 1: Tcj \fo\ versus Nc. in (161 (upper curve) compared to 
the numerical fit to the lattice results (lower curve) from [^ . 


p ~ p<d(9) A 2dep and pv k,- dtp (22) 


Inserting (22) into S yields 


5. Dyson Coulomb gas. We note that © coin¬ 
cides with the saddle point equation to 0 by re-writing 
it using Dyson charged particle analogy on with the 
energy 2S[z\ = J2i<j the pair interaction 


/ Jf\ 

- pI{S)W‘^Vp\) (23) 

with the potential 


G{zd = + oc{T)V{\zd) = G(0, - 9,) (17) 

At large Nc the ensemble described by 0 is sufficiently 
dense to allow the change in the measure. Following 
Dyson m we obtain 


W[p\ = 2a{T)[dep\s - TT!3{T)[dgp\H + dg 


/ dgp \ 

\pm) 


(24) 


For constant po and large Nc, (23) simplifies to 


Z[a, /3] —)■ 




with the effective action 


(18) 


m\,V2 / dtd9 {{dtpf - {dgpf 


(25) 


after the rescaling Vgt —)■ t with = 7rpo/3(r). (25) de¬ 


scribes sound waves in the large Nc space of holonomies. 


r[a,/3;p]= p{9)Gi9-9')pi9') 


8. Z(Nc) bubble. In a de-confined phase of infinite 
volume, the Yang-Mills ground state settles in one of the 
degenerate Z(Nc) vacua. In a finite volume, bubbles of 
different vacua may form |23j . Consider a de-confined 
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bubble of volume V2 immersed in a confined volume V2. 
In V2 all the eigenvalues are localized initially within a 
small A9 around the origin with p{t = 0,9) = N^jA9 = 
Pb, and zero otherwise. 

Using this piece-wise wave as an initial condition we 
solve (211 with A = 0 for simplicity. For large times r, 
the result is 


p{t,9) « Po - -pB sin 


A9 

2 


(26) 


which shows the relaxation of the piece-wise wave over a 
time r « l/vs set by the speed of sound. Using (26) in 


S yields the Euclidean action estimate for small A9 


Se{^ 2 ) 


V2 


(^moPB sin 



V2 (|lVemo)("27) 


The bubble formation probability or fugacity is e ®e(V 2) ^ 


9. Polyakov line in 1 -|- 3 dimensions. To extend 
our analysis to 1 -I- 3 dimensions, we approximate the 
Yang-Mills thermal state by a dense plasma of dyons and 
anti-dyons HIS]. This semi-classical description repro¬ 
duces a number of key features of the Yang-Mills phase 
both in the confined (center-symmetric) and de-confined 
(center-broken) phase. There are two key differences with 
the 1-1-2 dimensional partition function in (1). First 
the many-body energy 2S[z] = —21n4'o[z] in (1) is now 
shifted 




2S[z] ^ 2S[z] - 7 (T) [|(0,+i - 9i) -^ 


(28) 


with 7(r) = AttNc/Vs and / = the dyon fu¬ 

gacity m- Second and more importantly /3(T) = 2 and 
is not extensive with the spatial 3-volume V 3 . Finally, 
a{T) = T^V^/O. Since (6*^+1 — 9i) « l/2TTp{9i), then in 
the continuum the additional string of factors in (28) is 


Na 


J|(6i,+i -9i)^ ^ ^^fdep(e)Hi/27rp(e)) ^39) 


With this in mind, a re-run of the preceding arguments 
yields the Hamiltonian in (3]|4| with the shifted potential 


A ^ H + £^e-^«[^laelnp(0) (30) 

and A^cln7o[p] = / d9p{9)\n{p{9)/Nc)■ The hydro-static 
equation (|^ now reads 


^ttph{9) - 2A{9) = (^1 + j de\np{9) (31) 


The (3 = 2 contribution is now sub-leading and can be 
dropped. The corresponding solution to (31) is a local¬ 


ized density for ttci = d9p{9)cos9 ^ 0, and a uniform 
density po = Nc/2-k for ci = 0. Specifically 


pM. 

Po 


SttQTQ CO^d 

e 


(32) 


with c' = Ci/Nc and 7' = 'y/N((. The two parameters 77 = 
87ra(r)/7' and x = c' 7779 are fixed by the transcendental 
equations 


Il{x) TTX 
loix) 7770 


and 



27r^a; 

77 


(33) 


A solution exists only for 7' < 2a(T)l'K. Else the den¬ 
sity is uniform. Thus the transition temperature from 
center symmetric (confining) to center-broken (deconfin- 
ing) occurs for Q!(Tc)/ 7 (Tc) = l 2 .Nl '^i^h 

A = g^Ncl'3i'K^. For the dyon model, the fundamental 
string tension is given by a\ = (IVc/tt) sin(7r/iVc) A^/A [1]. 
Thus the model independent ratio in 1 -|- 3 dimensions 




37r 


v^Ti V8A^cSin^(7r/7Vc) 


(34) 


(34) compares favorably to the lattice result s iZ2\ even 
for small Nc as shown in Fig. At large Nc, (34) is con¬ 
sistent with the value of •y/3/27r in the string model [TB] . 



FIG. 2: Tc/v^oT versus in (34l. The dots are the lattice 
results from 1221. 


10. Conclusions. The hydrodynamical description of 
the Polyakov line captures aspects of the center dynamics 
in Yang-Mills theory in terms of the gauge invariant den¬ 
sity of eigenvalues. The hydro-static equations yield solu¬ 
tions that interpolate between a center symmetric (con¬ 
fining) and a center-broken (de-confining) phase. The 
transition temperatures normalized to the string tension 
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compare well to the lattice results over a broad range of 
Nc, and asymptote the string model results at Nc = oo. 
The hydrodynamical set-up supports a hydrodynamical 
instanton that describes the stochastic relaxation of the 
eigenvalues of the Polyakov line viewed as a fluid. The 
fluid supports sound waves that can be used to estimate 
the probability of formation of Z(Nc) bubbles. The re¬ 
laxation of a fluid of holonomies across the critical tem¬ 
perature may prove useful for understanding the onset of 


equilibration in a Yang-Mills plasma. 
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